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No Radio Afterglow from the Gamma-Ray Burst of February 28, 1997 

D. A. Frail 1 , S. R. Kulkarni 2 , D. S. Shepherd 2 , E. Waxman 3 

ABSTRACT 

We present radio observations of the gamma-ray burster GRB 970228 made with the 
Very Large Array (VLA) and the Owens Valley Radio Observatory (OVRO) spanning 
a range of postburst timescales from one to 300 days. A search for a time- variable radio 
source was conducted covering an area which included a fading X-ray source and an 
optical transient, both of which are thought to be the long wavelength counterparts to 
the gamma-ray burst. At the position of the optical transient sensitive limits between 
10 /iJy and 1 mJy can be placed on the absence of a radio counterpart to GRB 970228 
between 1.4 and 240 GHz. We apply a simple formulation of a fireball model which 
has been used with some success to reproduce the behavior of the optical and X-ray 
light curves. Using this model we conclude that the radio non-detections are consistent 
with the peak flux density of the afterglow lying between 20-40 fiJy and it requires 
that the optical flux peaked between 4 and 16 hours after the burst. 

Subject headings: gamma rays: bursts - radio continuum: general 



1 National Radio Astronomy Observatory, Socorro, NM, 87801, USA 

2 Division of Physics, Mathematics and Astronomy 105-24, Caltech, Pasadena CA 91125, USA 
institute for Advanced Study, Princeton, NJ 08540, USA 



- 2 - 



1. Introduction 

The gamma-ray burst of 28 February 1997 was a turning point in our understanding of these 
enigmatic objects, resulting in the first-ever discovery of X-ray and optical counterparts to a burst. 
Within the original 3-arcminute localization provided by the Wide Field Cameras (WFC) on board 
the BeppoSAX satellite, a previously unknown X-ray source 1SAXJ0501. 7+1146 was detected by 
the Narrow Field Instruments (NFI). Eight hours after the burst the flux of 1SAX J0501. 7+1146 
was 2.8xl0 -12 erg cm~ 2 s _1 (2-10 keV), but three days later its flux had dropped by a factor of 
20 (Costa et al. 1997). A comparison of V- and I-band images taken on 28 February and 8 March 
revealed an optical transient within a reduced error box, defined by the intersection of the WFC 
circle, the ±50" circle of the NFI, and the Interplanetary Network (IPN) annulus (van Paradijs et 
al. 1997, Hurley et al. 1997). 

Predictions of long-lived afterglows from gamma-ray bursts at X-ray, optical and radio 
wavelengths have been made for some time (e.g. Paczyhski & Rhoads 1993, Meszaros, Rees, & 
Papathanassiou 1994, Katz 1994, Meszaros & Rees 1997). In particular, a power-law decay in the 
observed long-wavelength flux with time is a generic consequence of a class of models known as 
"fireballs" . The gamma-rays are thought to be produced when the relativistically expanding blast 
wave (aka fireball) is slowed down by the ambient gas or is self-shocked by its own ejecta. The 
fireball accelerates particles in a shock which then radiate via the synchrotron process (Meszaros 
et al. 1994, Waxman 1997a, 1997b, Sari, Piran, & Narayan 1998). As the fireball expands it cools, 
shifting the peak in the spectrum to lower energies and resulting in delayed emission at longer 
wavelengths. 

The optical and X-ray decay from GRB 970228, as well as that from several other subsequent 
GRBs, is consistent at least to first order with one of the simpler formulations of these models 
(Meszaros & Rees 1997). Costa et al. (1997) fit a t 5 decay to the X-ray data with 5 ~ -1.33 + 0.12, 
while global fits to the X-ray, optical and infrared data give S = —1.2 (Wijers, Rees and 
Meszaros 1998) and 5 = —1.09 ± 0.23 (Reichart 1997). More recent optical fits, aided by a long 
time-baseline, yield 5 = -1.12 ± 0.08 (Garcia et al. 1997), 5 = -1.21 ± 0.02 (Masetti et al. 
(1997), and 5 = —1.10 ± 0.04 (Galama et al. 1998). The character of this decay agrees well 
with a fireball produced by a one-time impulsive injection of energy in which only the forward 
blast wave efficiently accelerates particles (i.e. the adiabatic piston model). The adiabatic model 
predicts a simple relation between the slope of the temporal decay 5 and the flux spectrum /?=|<5. 
Observationally the spectrum is not well determined, with [3 = —0.7 + 0.1 optically (van Paradijs 
et al. 1997) and X-ray values of j3 ~ —0.9 with large scatter (Frontera et al. 1998). The slope 
of the particle spectrum p (where (3 = (p — l)/2) inferred from these values of temporal decay is 
p ~ —2.6, not an unreasonable value for a relativistic shock (Blandford and Eichler 1987). 

Given the early success of the fireball model in predicting the gross properties of the optical 
and X-ray behavior of GRB 970228, it is reasonable to look for delayed radio emission from 
this burst. This is particularly relevant in the light of the discovery of the radio afterglow from 



-3- 



GRB 970508 (Frail et al. 1997a) which exhibited temporal and spectral behavior consistent with 
this model (Waxman, Kulkarni & Frail 1997). The properties of the fireball for GRB 970228 are 
well constrained by the optical and X-ray data. Thus the presence or absence of radio emission at 
late times has a bearing on the validity of this model. With this in mind, we began a radio search 
centered on the optical transient detected by van Paradijs et al. (1997). This Letter is a summary 
of our monitoring program for the first 300 days. 

2. Observations and Results 

Radio observations of a field centered on GRB 970228 began with the Very Large Array (VLA) 
on 1997 March 1.03 UT, 22 hours following the burst. Initially the VLA observations were made 
at a frequency of 1.43 GHz (A=21 cm), where the field of view is large enough to accommodate 
the preliminary 10' radius error circles from BeppoSAX. As the position of GRB 970228 became 
better known we increased the observing frequency first to 4.63 GHz (A=6.5 cm) and then to 8.46 
GHz (A=3.5 cm), beginning three days after the first VLA observation. On 1997 March 16.09 
UT a full synthesis single track was made with the millimeter interferometer at the Owens Valley 
Radio Observatory (OVRO) at 92.7 GHz (A=3.2 mm) and 242.2 GHz (A=1.2 mm). See Shepherd 
et al. (1998) for more details on the OVRO observation. 

The results of our VLA and OVRO observing campaign of GRB 970228, spanning postburst 
timescales from one to 300 days, are summarized in Table 1 and are plotted in Figure 1. The dates 
of each observation are given followed by the number of elapsed days, the observing frequency, the 
synthesized beam size and the rms noise on each image. The last column in Table 1 quantifies 
our principal scientific result, that between one and 300 days after the initial gamma-ray event no 
radio emission was detected from GRB 970228. One radio source was detected in the ±50" radius 
of the NFI (Frail et al. 1997b), but it lies outside the reduced IPN/WFC/NFI error box and has 
spectral and temporal characteristics that are consistent with it being an unrelated background 
extragalactic object. At the location of the van Paradijs et al. (1997) optical transient sensitive 
la limits between 10/xJy and 1 mJy can be placed on the absence of a radio counterpart to 
GRB 970228 between 1 and 250 GHz. The limits at 8.46 GHz are the most severe, with typical 2a 
limits of 100 /uJy. Adding the 8.46 GHz data from all 14 epochs together produces a 2a upper limit 
of 9 fiJy. The Westerbork Radio Synthesis Telescope (WSRT) and the Berkeley-Illinois-Maryland 
Array (BIMA) also observed GRB 970228 at 5 GHz and 86.4 GHz, respectively, producing upper 
limits a factor of two or more above those in Table 1 (Groot et al. 1997, Smith et al. 1997). 



3. Discussion 
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3.1. The Adiabatic Piston Model 

In this section we will briefly review the spectral and temporal behavior for the impulsive, 
forward shock dominated fireball model. It is treated in far more detail by Meszaros & Rees 
(1997), Waxman (1997a, 1997b) and Sari et al. (1997). We will not consider radiative models 
nor treat the case of rapid cooling (Vietri 1997b). The optical light curve of GRB 970228 has 
maintained its power law slope with no break even six months after the initial detection (e.g. 
Garcia et al. 1997), suggesting that the adiabatic model is an appropriate description of this 
system. Moreover, the absence of such a break suggests that the fireball is not far from being 
spherically symmetric (i.e. that it is not a jet). 

Initially the fireball ejecta have a bulk Lorentz factor r=10 2 -10 3 which falls as t -3 / 8 , where 
t is elapsed time in the observer reference frame. The blast wave ahead of the ejecta accelerates 
a power-law distribution of electrons with a minimum Lorentz factor 7 m and a power law index 
of p. The emergent spectrum in the co-moving frame has a synchrotron break at a frequency 
v m cc (where B' is the co-moving magnetic field) . Most of the radiation is emitted near the 

synchrotron break frequency. Above v m the spectral slope /3=(1— p)/2 is set by the shock processes. 
Below v m the spectral slope a is predicted to be 1/3 from synchrotron theory (Pacholczyk 1970, 
Katz 1994) but it could be flatter if there are a range of physical conditions (i.e. B', T) in the 
shock. 

The delayed emission at X-ray, optical and radio wavelengths originates from the deceleration 
of the fireball, which produces a shift of v m to longer wavelengths. The X-ray and optical 
flux will reach a maximum F m when v m moves into their respective wavelength bands with a 
time-dependence v m cxt _3//2 . The flux when v < u m rises as t 3a / 2 to a peak of F m , thereafter 
decaying as t 3 ^/ 2 when v > u m . There are two significant points to note here, first the frequency 
v m is independent of the fraction of fireball energy initially carried by rest mass and of the density 
of the medium into which it is expanding, and second, the peak flux density F m is the same at all 
frequencies. We further note, that as a general result it is easy to show that the radio emission 
increases until T ~a few, when the shock is only weakly relativistic. Since all the terms in u m 
depend linearly on T, u m cx T . At the time of the gamma-ray burst T ~ 300 and u m ~ 10 20 Hz, 
so scaling to the radio band [y m ~ 10 10 Hz) we obtain r ~ 1. 

The above behavior assumes that the emission emerges from an optically thin fireball. 
Synchrotron self absorption was predicted (Katz 1994, Katz & Piran 1997) to be important at low 
radio frequencies during the first few weeks of the fireball's existence. For GRB 970508 Shepherd 
et al. (1998) observed a low frequency turnover radio spectrum of the afterglow during the first 
week, consistent with synchrotron self absorption with an optical depth of one near 1.6 GHz. 
Frail et al. (1997a) observed a gradual rise in the 1.43 GHz flux at later times consistent with 
the fireball becoming optically thin at this frequency. The spectral slope in the self-absorbed part 
of the spectrum is v 2 when the shock is relativistic but becomes the familiar v^l 2 when T ~ 1 
(Katz 1994, Waxman 1997b, Katz & Piran 1997). Early estimates erroneously used a synchrotron 
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self-absorption frequency va of ^10 13 Hz (e.g. Meszaros & Rees 1997) but later corrections have 
shifted va to between 1-5 GHz (Waxman 1997b). In this paper it is pointed out that the addition 
a low energy population of electrons below 7 m can increase the synchrotron opacity without 
increasing F m , provided that these electrons are small in number compared to the high energy 
population. While this increases va somewhat, the fact remains that beyond 1-2 weeks after the 
burst the fireball emission should be optically thin at 8.46 GHz where most of our observations 
have been made. 



3.2. The radio data and constraints on the fireball model 

To first order the model described in §3.1 gives a fairly good description of the behavior of 
the X-ray and optical afterglow from GRB 970228. We will now use it to predict the late-time 
behavior of the afterglow at radio wavelengths. A simple-minded extrapolation of the spectral 
slope a, measured at the time when the optical and X-ray transients were detected, would predict 
flux densities in the centimeter range of 10-100 mJy Such flux densities were clearly not seen so 
we can infer that at this time the radio band v r was below v m . Frequencies below 5 GHz may 
have been initially strongly self absorbed (i.e. v r < v a i> < v m ) but this would not be the case for 
the majority of our data at 8.46 GHz nor our OVRO data at millimeter wavelengths (Waxman 
1997b). 

The radio emission peaks at F m at a time t r , delayed from the optical peak given by 
t r =t (v r /v o y 2 / 3 . For ty=8.46 GHz t r ~ 1400 t c for v r =92.7 GHz and t r ~ 300 t G . Neither F m 
nor t G are well constrained by the X-ray and optical afterglows since a distinct turnover was not 
observed. One approach used by Wijers et al. (1997) is to use the properties of the gamma-ray 
burst to constrain the values of F m and t G for the fireball. There is no a priori reason to assume 
that the GRB and the fireball are directly related in this way. In fact, it is currently held that the 
two are physically decoupled from one other (e.g. Meszaros 1997). 

Nevertheless, if we adopt t G ~ 670 s for R band used by Wijers et al. we find that the 
predicted radio emission from GRB 970228 would peak at ~ 10 6 s, thereafter decaying in the same 
manner as the optical and X-ray. At 8.46 GHz we have sensitive upper limits for the absence 
of radio emission at the location of the optical transient, spanning approximately two orders 
of magnitude in time (10 5 -10 7 s). For the value of F m ~ 5000 //Jy and the t _L2 decay fit for 
v > v m from Wijers et al. (1997), the VLA would have been expected to detect a bright radio 
source during its entire monitoring interval (>100 /uJy). Our millimeter observations at OVRO 
are less constraining. At the time of these observations, nearly two weeks after the burst, the 
millimeter emission from GRB 970228 is expected to have decayed below detectablity and thus it 
cannot provide any useful constraints on the fireball model. Wijers et al. (1997) assume a = for 
GRB 970228, so the flux is constant and equal to F m until v = v m . A more conservative approach 
is to set a to the asymptotic value of 1/3, so that the flux rises as t 1 / 2 . This is still at odds with 
our radio data and we conclude that the estimates of F m and t used by Wijers et al. are incorrect. 
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The early optical light curve for GRB 970228 provides additional constraints on F m and t c . 
In the adiabatic model the optical flux at any time t provided v Q > v m (i.e. the decay phase) is 
given by F =F m (u /u m ) /3 , or F =F m (t/t ) 3//2/3 . At the same time provided v r < u m (i.e. rising 
phase) then F r =F m (v r /v m ) a or F r oct 1 / 2 for a = 1/3 (Waxman 1997a, 1997b). We can solve these 
two equations to relate the rising radio flux density at any time in terms of the decaying optical 
flux at the same time. 

F r = F (u r /u ) a (t/tj 3 / 2 ^ 

As noted earlier the radio flux at 8.46 GHz is expected to rise to a maximum F m at a time 
t r ~ 1400 t thereafter it would decay identically to the optical and X-ray. Small values of t give 
rise to large values of F m . 

A compilation of flux measurements for the optical transient are given by Galama et al. 
(1997) and updated in Galama et al (1998). The earliest measurement is by Pedichini et al. (1997) 
on February 28.81 UT which was translates to an R-band magnitude of 20.5±0.5 (Galama et al. 
1997). A similar but more accurate value of 20.9±0.14 was measured on February 28.99 UT by 
van Paradijs et al. (1997). Converting these R-band magnitudes to flux densities we derive values 
of 18+ 41 [iJy and 12±2 fi3y, respectively. No clear signature of a turnover has been seen in the 
optical data and in fact the power law decay of GRB 970228 has maintained its character with 
small deviations beginning 16 hrs after the burst to the present time (Galama et al. 1998). Thus 
the peak flux density must have been prior to this time, for which we derive F m >18 fiJy and 
t < 6 x 10 4 s. 

The radio upper limits in Table 1 can constrain F m and t Q further still. Using the fit to 
the optical light curve from Galama et al. (1998) and the above equation, we plot a number of 
hypothetical radio light curves in figure 1 for t o =6xl0 2 s, 6xl0 3 s, and 6xl0 4 s. The absence 
of readily detectable radio emission 300 days after the original event requires t D >2.8xl0 4 s and 
F m < 40 /LtJy. Similarly, an upper limit on t D is set by the date of the earliest detection of optical 
emission from GRB 970228 of 18 ^Jy at 6x 10 4 s by Pedichini et al. (1997). If the adiabatic fireball 
model is to remain a valid description of the afterglow emission from GRB 970228 then it predicts 
that the radio afterglow will peak around 40 to 18 /txJy between 2xl0 7 s and 8xl0 7 s. Castander 
& Lamb (1998) have argued in favor of visual extinction Ay=1.3 towards GRB 970228, a value 
some two to three times larger than used in most papers to date. This has the effect of raising 
F m , making the radio non-detections even more problematic. The radio emission may be weaker 
than these predictions when the decelerating fireball eventually becomes sub-relativistic (Wijers 
et al. 1997). However, such a transition should first be seen in the optical light curve as a break 
in the power law slope and no such break has been seen thus far. 
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4. Conclusions 

The detection of an optical transient, thought to be the afterglow from GRB 970228, has made 
it possible to perform a search for time-variable radio emission at the optical position. Delayed 
emission at longer wavelengths (X-ray, optical and radio) is a generic prediction of all fireball 
models. VLA and OVRO observations have been made that span a range of postburst timescales 
from one to 300 days putting limits on the absence of a radio counterpart to GRB 970228 between 
10 fj,Jy and 1 mJy Applying a simple version of the fireball model which has been used successfully 
to fit the temporal behavior of the decaying X-ray and optical emission from GRB 970228 suggests 
that the radio afterglow has yet to peak. If this is correct then continued deep imaging in the 
coming months offers the possibility for the detection of a weak but increasing radio signal. 
Detecting this emission would be a powerful verification of the fireball model as described by 
Meszaros & Rees (1997). 

Indeed, as the recent detection of GRB 970508 has taught us, radio afterglows yield unique 
GRB diagnostics that are not obtainable by any other means (Frail et al. 1997b). Unlike optical 
or X-ray wavelengths one is presented with the possibility of following the full evolution of the 
fireball emission through all its different stages; first while it is optically thick, then as it slowly 
rises to a peak flux density and thereafter decays, making a transition from an ultra-relativistic to 
sub-relativistic shock. Furthermore, both the scintillation of the radio source (Goodman 1997) 
and its flux density, when it is synchrotron self-absorbed (Katz 1994), allow a determination of the 
size and expansion of the fireball. If the radio afterglow is going to be detected from GRB 970228 
and others like it, then continued long-term monitoring is going to be required at the microJansky 
level. 

The VLA is a facility of the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. SRK's research is supported by the National Science 
Foundation and NASA. DAF thanks Ralph Wijers for useful discussions. 
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Figure 1.— Upper limits on the flux density of a radio counterpart to the optical transient 
associated with GRB 970228. VLA points are given by the filled triangles (1.4 GHz), a cross 
(4.6 GHz) and the filled squares (8.5 GHz). OVRO points are the filled circles at 240 GHz (top) 
and 90 GHz (bottom). The dashed lines are the expected radio light curves at 8.46 GHz for the 
adiabatic piston model using the observed optical decay (Galama et al. 1998) scaled to the radio 
and assuming several values for the timescale of the onset of the decay phase (see text for more 
details). 
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